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On Cosmic Censorship in String Theory 
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We investigate the possibility of cosmic censorship violation in string theory using a characteristic 
double-null code, which penetrates horizons and is capable of resolving the spacetime all the way 
to the singularity. We perform high-resolution numerical simulations of the evolution of negative 
mass initial scalar field profiles, which were argued to provide a counterexample to cosmic censor- 
ship conjecture for AdS-asymptotic spacetimes in five-dimensional supergravity. In no instances 
formation of naked singularity is seen. Instead, numerical evidence indicates that black holes form 
in the collapse. Our results are consistent with earlier numerical studies, and explicitly show where 
the 'no black hole' argument breaks. 



I. INTRODUCTION 

Cosmic censorship conjecture remains an unsolved 
problem in general relativity. In the absence of proof, 
finding an acceptable counterexample is very important, 
as it would resolve the issue one way or the other. But 
even that remains elusive. While there are numerous ex- 
amples of initial conditions that form a naked singularity 
in general relativity, none of them are generic as required 
by the terms of the cosmic censorship conjecture. 

Recently, Hertog, Horowitz, and Maeda argued that 
cosmic censorship is generically violated in asymptoti- 
cally anti de Sitter spacetimes in particularly in five- 
dimensional supergravity arising from string theory con- 
structs The counterexample consists of an initial 
data set which is known to form a singularity and an 
argument that it cannot form a black hole in the col- 
lapse. That would leave the possibilities that singularity 
is naked, which would disprove cosmic censorship con- 
jecture, or consumes the entire spacetime (big crunch), 
which initially was deemed unlikely by the authors. 

These papers have caused quite a controversy in the 
literature. The arguments of were called into ques- 
tion by a number of papers, on grounds of both analytical 
0, 1^] and numerical ^, 0, calculations. The authors 
themselves pointed out a gap Q in their first example 
, and reconsidered the possibility of a big crunch Q . 

The present paper investigates the possibility of cos- 
mic censorship violation numerically. Earlier numerical 
studies IHSH suggestive of a black hole formation. 
However, they are not conclusive due to the flawed coor- 
dinate choice, which breaks at the moment of apparent 
horizon formation, and prevents the code from seeing the 
whole spacetime, thus leaving the possibility of a singu- 
larity subsequently becoming naked. This shortcoming 
is resolved in the present paper by using a characteristic 
double-null code. 
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Conceptually, the code is similar to the one used to 
investigate the global structure of the spacetime in the 
critical collapse of the scalar field fl^ , although the im- 
plementation is entirely new. Several issues specific to 
asymptotically AdS spacetimes had to be addressed, in 
particular boundary conditions at infinity. 

We focus on cosmic censorship counterexample in the 
string theory context Q , and perform high-resolution nu- 
merical simulations of the evolution of two families of 
negative mass initial scalar field profiles for various val- 
ues of parameters and different boundary conditions at 
the cutoff. In no instances formation of naked singularity 
is seen. Instead, numerical evidence indicates that black 
holes form in the collapse, in agreement with Our 
results clearly show where the 'no black hole' argument 
of 2] fails. 

The paper is organized as follows: In Section^ the 
issue of a coordinate choice is discussed in view of the 
requirements of testing cosmic censorship conjecture nu- 
merically. The scalar field evolution equations and the 
implementation of the algorithm are described in Sec- 
tion The numerical results of investigation of cosmic 
censorship in string theory are given in Section llVI and 
summarized in Section Ivl 



II. COORDINATE CHOICE 

In numerical investigations of singularity formation 
and global structure of the spacetime, the coordinate 
choice plays a crucial role. To study the global struc- 
ture, coordinates must cover the entire spacetime. To 
see the singularity inside a black hole, coordinates must 
penetrate the horizon. To prevent the singularity from 
corrupting the rest of the spacetime once it forms (if it 
is not naked, that is), the information propagation speed 
(both numerical and true characteristic) must be under 
control. These issues cannot be avoided if one wishes to 
study cosmic censorship numerically. 

In this paper, we adopt double-null coordinate system, 
which is particularly suited for the above requirements. 
The spherically symmetric (n-|-2)-dimensional spacetime 



metric can be written as 



(1) 



where dJl^ is the metric of a unit ri-climensional sphere, 
and the two-manifold metric 1^7^ = e^^^diy^ is written 
in exphcitly conformally-flat form. Throughout the pa- 
per, we will be freely switching between Minkowski and 
double-null coordinates on the flat two-manifold 



dr]^ 



-di^ + dx^ 



-A du dv. 



(2) 



with u = [t — x)/2 and v — {t + x)/2. The familiar form 
of anti de Sitter metric in spherically symmetric static 
coordinates 



ds' = - 1 



dt^ 



dr"^ 



'dni 



transforms to 



r = I tan — , cr = in cos — 



(3) 



(4) 



in coordinate system . One should note that the pecu- 
liar property of AdS that the spatial infinity "is not very 
far away" is reflected in the fact that a finite interval 
X e [0, |- ^] covers the whole AdS spacetime in conformal 
coordinates (0. Indeed, these are the coordinates used 
in construction of a Carter-Penrose conformal diagram of 
the AdS spacetime. 

In spherical symmetry, one can define a local mass by 
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2m 



(5) 



The function / carries information about the spacetime 
structure: it is negative in the trapped (or anti-trapped) 
region, positive in the regular region, and vanishes on the 
apparent horizon. The mass function defined by |SJ| coin- 
cides with ADM mass in asymptotically flat spacetimes. 
In asymptotically AdS spacetimes, however, it diverges 
near spatial infinity. The physical cause for this is simple: 
AdS has a constant (negative) energy density while the 
volume is infinite. To avoid infinities, one can subtract 
the divergent part from the definition of the mass 



2 j,n+l 



(6) 



where £ is the curvature radius of the asymptotic AdS. 
We will refer to this mass definition as the reduced mass. 

Although traditionally used in numerical relativity, 
Schwarzschild coordinates, like the ones in metric (O, are 
not suited to study the possible violation of cosmic cen- 
sorship in dynamical evolution. They do not penetrate 



horizons (at the apparent horizon, / 



9' 



— 0, so the 



metric coefficient grr necessarily diverges), and therefore 
have no chance of seeing, for example, the destruction of a 
black hole by infall of a negative mass and "baring" of the 
singularity inside (if it actually happens) . Unfortunately, 
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FIG. 1: Coordinate choice. 
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Schwarzschild coordinates were the coordinates used in 
recent numerical studies 0, 0, of the possible coun- 
terexamples to cosmic censorship in asymptotically AdS 
spacetimes 0, Q • Although the formation of a trapped 
surface is seen, the further fate of the spacetime remains 
unknown. This renders the results of 5, 6, 7] largely in- 
conclusive as far as the possible violation of the cosmic 
censorship goes. 

There are coordinate systems other than which 
do penetrate the horizon, most notably the ingoing 
Eddington-Finkelstein coordinates 



ds^ = -fe^s^v^ + 2eBdvdr + r'^dflf 



(7) 



However, in implementation of the evolution code based 
on this coordinate system (which we did as well, although 
it is not going to be further discussed in this paper) 
one encounters a problem that the numerical information 
propagation speed is higher than the true speed of the 
outgoing characteristic (which varies on the grid). This 
leads to termination of the numerical evolution in the re- 
gions causally disconnected from the singularity, which 
is not a satisfactory behavior for a code that aspires to 
resolve the global structure of the spacetime. 

The true strength of double-null coordinate system 
is the advantage of knowing the speed of information 
propagation everywhere in advance. The importance of 
this point cannot be overemphasized, as numerical infor- 
mation propagation properties is what makes or breaks 
the numerical code which has to deal with singularities 
before any other factor even comes into play. 

The final issue that needs to be discussed before we 
go over to evolution equations is the gauge choice. The 
form of the metric (^3) is left invariant under rescaling of 
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null coordinates 



U{u), 



V{v). 



(8) 



We use this gauge freedom to place initial and boundary 
surfaces at known coordinate locations, namely fix center 
of spacetime and constant-r cut-off at a: = and x — c 
correspondingly, and put initial data surface at t = 
timeslice, as illustrated in Fig. ^ This can always be 
achieved simply by putting u = w at r = and u = —v 
at initial surface, while rescaling v to have v — u — c at 
cut-off. In fact, it does not even fix the gauge uniquely. 
The residual gauge freedom © is determined by 

V{u) ~ U{u) = 0, u > 0, (9a) 
V{u + c) - U{u) = c, u> -c, (9b) 
V{-u) + U{u) = 0, > M > -c, (9c) 

and amounts to an odd periodic function lo of period c, 
with V{u) — U{u) = u + Lo{u). 



are recovered by variation with respect to the (flat) two- 
metric "qab- The traceless part of a symmetric two-tensor 
has two independent components, which means there are 
two constraints. They can be taken either as {uu\ and 
{vv} or {xt} and {xx} -\- {tt} components of ((T^ . 

The evolution code uses a semi-constrained algorithm. 
Rather than solve dynamical equation ()12cfl directly, in- 
going [uu) null constraint equation 



2r, 







(14) 



is integrated to obtain one of the metric coefficients. The 
integration is particularly easy to do if one defines a new 
field variable 

-2a-ln(-r„), (15) 
in terms of which the constraint H14(l is written simply as 



r < 

9,u = -- 



(16) 



III. SCALAR FIELD EVOLUTION 

The dynamics of a scalar field with potential V{(j)) 
minimally coupled to gravity in iV = 2 + n dimensions is 
described by the action 
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{i^-g^-- -2y(0)}y3^d^^. (10) 



Spherical symmetry effectively reduces the problem to 
evolution of three non-linearly coupled scalar fields on a 
flat two-dimensional manifold. Substituting the spher- 
ically symmetric metric into the scalar field action 
(|10|l and integrating over the n-dimensional spherical 
subspace, one obtains the reduced action 



^ oc ||li(^[(Vr)2+e-2-] -^(Vr-Vfi) 

-(V0)2-2e"2^1/(0)}r"d23;^ (11) 

which describes the dynamics of a spherically symmetric 
gravitating scalar field. Here and later the differential 
operators (gradient V and D'Alembertian □) are taken 
with respect to the flat two-dimensional metric (O . Vari- 
ation of the reduced action with respect to field variables 
a, 4>, and r gives equations of motion 



l + i^rY = (1-41^ j, (12a) 

(12b) 



□0 -I- - (Vr • V0) = e^^'^F', 
2 r 2 n 



(12c) 



while constraint equations 

Traceless [r^ab + 2 r(^<jO-,b) + ^ 0,a(/',b] = (13) 



The field variable g actually has a meaning beyond a mere 
computational convenience. It is the same g which enters 
coefficients of the Eddington-Finkelstein metric (TJ. 

Equation H12a|) warrants further discussion. It involves 
a delicate cancellation of terms in both small r and large 
r regimes, which makes it susceptible to discretization 
errors. The loss of precision can be catastrophic, leading 
to numerical instability developing at either r = or 
large r and subsequent code failure. After trying various 
discretization schemes, we settled on discretizing a form 
of equation H12a|) which eliminates non-linear gradient 
term in favour of a box operator for small r 



□ (r 



n{n — 1) 



1 - 



2r^V 
n(n— 1) 



(17) 



while for large r, the equation H12al) re- written in terms 
of inverse radius p — 1/r is discretized instead 

pUp - (n + l){Vpf = -p,.e3 [{n - l)p^ - ^V] . (18) 

This largely circumvents the stability problems men- 
tioned, while still being possible to discretize efhciently 
The Cauchy problem for evolution equations (|12l) re- 
quires specification of six functions of one variable at the 
initial spacelike surface: the values of three fields r, cr, </> 
and their time derivatives r, ct, (j). These six functions are 
not independent in general relativity; they must satisfy 
two constraint equations (|13|l . We restrict our attention 
to time-symmetric initial data. In covariant form, the 
requirement of time symmetry is written as 



Kab - 0, n • V0 = 0, 



(19) 



where K-ab is an extrinsic curvature and n is the normal to 
the initial data surface. In our gauge, the initial surface 
is t = 0, so the condition of time symmetry is simply 



4 



that the time derivatives f, ct, (j) vanish. This requirement 
satisfies the {tx} constraint identically, which leaves three 
functions subject to one constraint for specification of the 
initial data. Of the two freely specifiable functions, one 
is physical - the initial scalar field profile, and the other 
is a gauge choice fixing residual gauge freedom 0. 

Once the scalar field profile is specified on the initial 
time slice, the rest of the variables are obtained by inte- 
grating 



2m . „ 

1 1 p3 



5,r 



(20a) 

(20b) 
(20c) 



Equation H20a() follows from definitions (|5I15() . Equation 
(|20b|l is the remaining constraint equation rc-written in 
terms of the physically interesting mass function m. Our 
gauge choice is implicitly given by equation (|20c() ; it cor- 
responds to the requirement g = on the initial time 
slice. The integration is performed numerically using 
fourth order Runge-Kutta algorithm with constant step 
size. 

Since in asymptotically AdS spacetimes the light rays 
take finite time to reach spatial infinity and reflect back, 
the issue of boundary conditions at infinity is a physical 
one, and cannot be avoided. Various boundary condi- 
tions for scalar field </) at infinity have been suggested, 
specifics of which we will discuss later. Additional com- 
plication for numerical evolution is that it is unfeasible 
to include infinity on the grid, so the spacetime has to 
be cut-off, which we do at a constant (large) radius rc, 
and impose boundary conditions for the scalar field there. 
The remaining boundary condition at the cut-off is the 
one for the metric coefficient 



5,t 



P.xt 
P,x 



2 p( 

n 



(21) 



which follows from the {tx} constraint equation H13() . In 
the center of the spacetime, one has the usual regularity 
condition for the scalar field 



Vr • = 0. 



(22) 



The numerical evolution scheme is illustrated in Fig. |21 
The field variables are discretized on a square spacetime 
grid of spacing e. Discretization of the evolution equa- 
tions is done using the leapfrog scheme, which is second 
order accurate and non-dissipative. In particular, the 
gradient and box differential operators are discretized as 




FIG. 2: Numerical evolution scheme. 

Neumann boundary conditions (|22|) require asymmetric 
discretization of the derivatives for second order accuracy 



^'(0) = ^(4Xi - 3Xo - X2). (24) 

The code was implemented in Fortran, and is efficient 
enough to run on personal workstation class computers. 
Since only three subsequent time steps are stored, the 
code is not memory-limited even for very large grids. All 
the numerical results presented in this paper were ob- 
tained in simulations with the grid size of 65539 points. 



IV. COSMIC CENSORSHIP IN STRING 
THEORY 

Hertog, Horowitz, and Maeda argued that cosmic cen- 
sorship is generically violated in string theory In 
their counterexample, they considered A/" = 8 gauged su- 
pergravity in five dimensions, which is thought to be a 
consistent truncation of ten dimensional type IIB super- 
gravity on . They picked a single scalar component 
which does not act as a source for any of the other fields, 
thereby reducing the problem to the one of Einstein grav- 
ity in five dimensions minimally coupled to a scalar field 
(j) with negative potential 



-2e 



20/\/3 



4e 



(25) 



The potential is unbounded from below, with a maximum 
value of 1^(0) = —6 which corresponds to the asymptot- 
ically AdS spacetime with curvature scale £ = \. They 
constructed time-symmetric initial data which has nega- 
tive reduced mass /iq of the initial configuration, namely 



A 
A 



r<Ro 
r>Ro 



(26) 



A"rt — Alt 

2'e ' 



V(A 



Aup — Adn 

21 ' 



OX = ^(Aj-t -|- Alt — Aup — Adn)- 



They further argued that a black hole cannot form in the 
(23a) collapse of this initial configuration (as Schwarzschild- 
AdS black holes have positive mass), whereas it is possi- 
ble to show that singularity does form in the collapse of 
(23b) homogeneous part of the field, and hence concluded 
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that singularity must be naked and cosmic censorship 
conjecture is violated. 

We use the numerical code described in the previous 
section to solve Einstein and scalar field equations for 
the spherically symmetric collapse of the profile (|26|l . and 
determine the global structure of the resulting spacetime. 
We find that black holes do form in the collapse of initial 
data 1)26(1 , contrary to the argument of Hertog et. at. Q . 

We place the spacetime cut-ofF at Tc = 25, which is 
large enough so that less than one percent of the initial 
profile mass /io(oo) lies beyond the cut-ofF in all cases. 
At the cut-ofF, we impose either Dirichlet or Neumann 
boundary conditions for the scalar field 



r—^ — const, 

= T, = 0, 



[Dirichlet BC] 
[Neumann BC] 



(27a) 
(27b) 



Neumann boundary conditions (|27b|l correspond to the 
"standard" boundary conditions on the fall-off of the 
scalar field (/) = a{t)/r'^ in AdS-CFT correspondence, but 
are written in this form for computational convenience. 
We also try alternative suggestion for the initial field 
profile. 



, Inr 
A— , 



r <i?o 
r>Ro 



(28) 



which has a slower fall-off rate and much larger negative 
masses than (j^ . 

Typical spacetime obtained in numerical simulations 
is shown in Fig. |3| It corresponds to evolution of ini- 
tial field profile (|26|l with values of parameters A — 16, 
i?o = 4 (00 = 1-0) subject to Dirichlet boundary condi- 
tions H27a|l . The four panels show density plots of metric 
coefficients r{x,t) (top left) and ■m{x,t) (bottom left), as 
well as the scalar field (f>{x, t) (top right) and its gradi- 
ent g^'^ 4) ^^4) ,^ (bottom right). Contours on r{x,t) plot 
show isolines of constant radius. The upper thick line 
on m(a;, t) plot is the singularity (where spacetime evo- 
lution terminates), while the lower thick line is the ap- 
parent horizon (found by condition / = 0). Due to our 
coordinate choice Q , the plots in Fig. O essentially are 
Carter-Penrose conformal diagrams, and show the global 
structure of the spacetime. It is clear that the black hole 
has formed in the collapse. 

To see if the black hole settles into an almost static 
configuration or continues to evolve, we plot the size of 
the apparent horizon tbh as a function of advanced time 
V in the left panel of Fig. 0] After the stage of growth 
corresponding to collapse of the homogeneous part of the 
initial scalar field profile, and sharp step-like increase in 
size when the outgoing portion of the wavepacket reflects 
from infinity and falls in (these features in the scalar 
field profile are most clearly seen in gradient plot in the 
bottom right of Fig.O, the black hole settles to almost 
constant size. The apparent horizon (shown in bottom 



left panel of Fig.|21l approaches null direction towards the 
end of the evolution as well, which is consistent with the 
black hole settling into a static configuration after the 
formation. 

The formation of the black hole does not contradict 
the negative initial reduced mass fj, of the profile. For 
Dirichlet boundary conditions (|27a|l . the reduced mass 
at the cut-off is conserved, but the black hole has scalar 
hair which hides the positive black hole mass and makes 
the total mass at the cut-off negative. No-hair theorems 
[Til do not apply in this case as the cut-off is at a 
finite distance, while the proof of the no-hair theorems 
involves exclusion of the growing mode at infinity. For 
Neumann boundary conditions l|27bl) . black holes do not 
have scalar hair, but the reduced mass at the cut-off is 
not conserved, and grows positive to accommodate the 
formation of the black hole, as shown in right panel of 
Fig.H 

In search of violation of cosmic censorship conjecture, 
we have performed many high-resolution numerical sim- 
ulations for different initial profiles (|26|l and (|28|l . dif- 
ferent boundary conditions (|27a|l and l|27b|) . and various 
values of the parameters (j)o and Rq, but all in vain. In 
no cases formation of naked singularity was seen. The 
global structure of the spacetime remains similar to that 
of Fig. 13 with black holes of various sizes formed in the 
collapse. Our numerical results are summarized in Ta- 
ble m which shows sizes and masses of the black holes 
produced. In addition, the spacetime mass at the cut-off 
after black hole formation is shown for runs with Neu- 
mann boundary conditions. The quoted values should 
be considered approximate, as the number of useful grid 
points shrinks near the end of the evolution and precision 
suffers somewhat (especially for smaller black holes) . It is 
worth noting that the black holes formed in the runs with 
Neumann boundary conditions indeed have (almost) no 
hair, as the masses evaluated at the horizon and at the 
cut-off are almost the same. 

As Fig. |31 shows, the size of a black hole formed in 
the collapse of the truncated 1/r^ field profile (|26|l scales 

like tbh = cxRo'Po^'^ ~ otA^^"^ as a function of profile pa- 
rameters, with coefficient a « 0.627 for Dirichlet and 
a ~ 0.693 for Neumann boundary conditions. This 
agrees well with the lower-bound estimate based on an 
analysis of the collapse of the homogeneous part of the 
profile Although the latter underestimates the value 
of numerical coefficient (yielding a = 0.58), it is only by 
about 10-20%. 



V. CONCLUSION 

We have developed a double-null characteristic code 
implementing iV-dimensional spherically symmetric evo- 
lution of a minimally coupled scalar field with potential in 
asymptotically AdS spacetimes, and used it to study the 
possibility of cosmic censorship violation in string theory. 
No instances of formation of naked singularity were seen 
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FIG. 4: The size of the apparent horizon (left) and the reduced mass fj, at the cut-off (right) for the same initial conditions as 
in Fig. 121 Solid and dashed lines show the values corresponding to Dirichlet and Neumann boundary conditions at the cut-off. 
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truncated field profile 
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68.0 
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4.0 
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11.542 
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1.5 
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4.0 


-945.3 
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24.8 



TABLE L Summary of tfie numerical results from evolution of the truncated 1/r^ (left) and ln(r)/r'^ (right) initial field profiles 
given by equations 12611 and 1281 for various values of parameters. Black hole sizes and masses are approximate. 
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FIG. 5: The size of the black hole formed in the collapse of the truncated 1/r^ field profile (1261 as a function of parameters, 
shown for Dirichlet (left) and Neumann (right) boundary conditions at the cut-off. 



in high-resolution numerical simulations of the evolution 
of two families of negative mass initial scalar field profiles 
for various values of parameters and different boundary 
conditions at the cutoff. 

Our results indicate that black holes form and reach 
steady-state configuration in the collapse, despite the 
negative initial mass of initial data. Either the spacetime 
mass becomes positive (for Neumann boundary condi- 
tions, where mass /i is not conserved), or the black hole 
covers itself with negative mass hair (for Dirichlet bound- 
ary conditions, where no-hair theorems do not applybe- 
cause of the finite cutoff), as was pointed out in 

The possibility of naked singularity formation cannot 
be ruled out completely by results of numerical simula- 
tions, as they only sample a finite number of initial con- 
figurations, but in view of the above it seems unlikely. Al- 



though no evolutions leading to big crunch were observed 
from initial configurations considered in this paper, that 
possibility should be explored further in the future work. 

To summarize, the cosmic censorship conjecture was 
found to hold in all the examples we have studied. 
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